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In order to study the effects of strong magnetic field on Kondo insulators, we calculate magnetization 
curves and single-particle excitation spectra of the periodic Anderson model at half-filling under finite 
magnetic field by using the self-consistent second-order perturbation theory combined with the local 
approximation which becomes exact in the limit of infinite spatial dimensions. Without magnetic field, 
the system behaves as an insulator with an energy gap, describing the Kondo insulators. By applying 
magnetic field to f-electrons, we found that the energy gap closes and the first order transition from 
insulator to metal takes place at a critical field He- The magnetization curve shows a jump at H^. 
These are consistent with our previous study in terms of the exact diagonalization. Relationship to the 
experiments on YbBi2 and some other Kondo insulators is discussed. 
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1 Introduction 

The materials called as Kondo insulators show the 
Kondo behavior at high temperatures as typical 
for the heavy fermion metals, but develop energy 
gaps at low temperatures. Examples of such ma- 
terials are SmS,§ SmBg,!] YbBi2,§ TmSe,g 
Ce3Bi4Pt3,|P etc. The energy gap seems to be 
formed due to the hybridization of the 4f states 
with the conduction bands, but it is considered 
that the strong correlation is playing an impor- 
tant role in determining the size of the gap and 
the other properties|6|. FeSi is also considered to 
fall into this family[3|. CeNiSn[0 is believed to 
have a pseudogap, or even to be a semimetal at 
lowest temperatures. Q 

In order to know the characters of the energy 
gap of the Kondo insulators, a stimulating ex- 
periment was performed on YbBi2 by Sugiyama, 
et al.^ In their experiment, the application of a 
strong magnetic field of about 47T destroys the 
gap and the system undergoes a transition from 
insulator to metal. The magnetization starts to 
steeply increase at a critical field. 

It is an interesting theoretical issue to investi- 
gate how the strongly renormalized ground state 
of the Kondo insulators will be modified and how 
the gap is destroyed by the magnetic field. For 
these purposes, we previously applied the recently 
developed techniques to treat the strongly cor- 
related electron systems in the infinite spatial 
dimensions [p^ pT| . In this scheme, the study 
of the strongly correlated lattice system is re- 
duced to solving the impurity problem in an ef- 
fective medium self-consistently. Therefore, it 



is sometimes called as the local-impurity self- 
consistent approximation or the dynamical mean 
field theory. [0 We solved the effective impurity 
problem by the exact diagonalization (EXD) of 
the finite system for conduction electrons and 
calculated the single-particle excitation spectra 
to determine the energy gap.p^, p^ , Q (The 
paper 1 12 will be referred to as I in the following.) 
We found that the gap decreases linearly with the 
magnetic field and vanishes abruptly at a certain 
critical field He, which is of the order of the impu- 
rity Kondo temperature, and that the insulator- 
to-metal transition takes place. The magnetiza- 
tion curves exhibit a jump at He- These features 
suggest that the transition is of first order. 

We discussed in I that there are two effects 
of the magnetic field on the energy gap in the 
Kondo insulators. One is the Zeeman shift of 
the up- and down-spin bands, which leads to the 
closing of the gap. The other is the reduction 
of the strong renormalization due to correlation, 
which results in the increase of the renormaliza- 
tion factor z = (1 - dT.f{e)/de)-\=o (S/(e) de- 
notes the f-electron self-energy) from 2 ^ 1 to 
z ^ \. This gives rise to the increase of the gap 
since the gap Eg is estimated to be of the order 
of zV"^ /W , where V denotes the mixing matrix 
element and W the conduction band width. We 
found that at least at half-filling in the symmet- 
ric case z does not change by the field as far as 
h < he- Therefore, the energy gap closes merely 
by the Zeeman shift of the bands, although an 
abrupt transition to metal takes place before the 
gap completely closes in this way. This first order 
transition may be due to the complicated many- 



body effects. However, it might be possible that 
it would be due to the finite size effect of the con- 
duction electron system. 

Recently, Carruzzo and Yu||l^ and Tsutsui, et 
al.[l^ claimed that the gap remains finite even 
under strong magnetic field and the field-induced 
insulator-to-metal transition does not occur in 
one dimensional Kondo insulators. Their con- 
clusions still seem to be indefinite because of 
their numerical methods (density matrix renor- 
malization group and exact diagonalization, re- 
spectively). However, it is interesting to study 
whether such transition occurs or not in higher 
dimensions. 

In order to clarify these points, we will use in 
the present paper the self-consistent second-order 
perturbation theory (SCSOPT),0, |l| and 
reinvestigate the character of the field-induced 
insulator-to- metal transition in Kondo insulators. 
As will be explained below, we have found again 
that the transition does occur and is of first or- 
der, so that our previous conclusions mentioned 
above are correct at least qualitatively. For com- 
parison, a calculation by the iterative perturba- 
tion theoryl^ (IPX) was also performed, but we 
found that IPT yields quantitatively insufficient 
results under finite magnetic field. We will also 
comment on the possible relationship to the ex- 
periments on YbBi2 and some other Kondo insu- 
lators. 

2 Self-consistent Second-order Per- 
turbation Theory and Iterative 
Perturbation Theory 

We use the periodic Anderson model (PAM), 

ij(T iCT 

ia i 

where Suf^ = nfcr-{nf„), Ej^ = Ef+U{nf-^)- 
ah, h — fisH and H denotes the magnetic field. 
N is the number of lattice sites. The other nota- 
tions are standard. We neglect the orbital degen- 
eracy and assume the symmetric case Ef = —U/2 
together with the half-filling condition to express 
the simplest model to the Kondo insulators. Since 
the g-factors for f- and the conduction electrons 
{gc and 5/, respectively) are different, we neglect 



gc for simplicity and apply magnetic field only 
to f-electrons, although a more general treatment 
is possible. We also assume a paramagnetic 
ground state. Possible antiferromagnetic states 
are regarded as being suppressed by the intro- 
duction of an appropriate frustration. pl[ 
In SCS0PT[|7|, 

the self-energy is calculated as 
/(-ei)/(e2)/(e3) + /(ei)/(-e2)/(-e3) 



e — ei -I- £2 + £3 + «<5 

where Pa{f) = (— l/7r)ImG'/CT(e + i5) is the f- 
electron density of states, /(e) denotes the Fermi 
function and 



(2) 
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Here denotes the energy of conduction elec- 
trons. These equations are converted into the fol- 
lowing forms: 

Jo 

+ a^{t)P^a{t)a-a{t)], (4) 



acr(i) 



dee— V.(e)/(±e), (5) 



G/.(e) = g.(e) [l -I- V^g.{e)F{e - V'g.{e))] 



g,{e)^{e~Ef,~^f'y{^)r' 



(6) 
(7) 



and 



F{z) 



{2/W)[z/W~^{z/WY-l]. (8) 



The last equation was obtained for the semi- 
circular density of states of the conduction elec- 
trons: p(e) = {2/ttW)[1- .Jl-{e/WY\. These 
equations should be calculated self-consistently. 

In IPT,[0 the lattice problem is converted 
into solving an effective impurity problem, so 
that /Ocr(e) in eqs.(^ and (||) should be replaced 
with that for the local Green's function, Pa{^) — 
(— l/7r)ImG/cr(e + i5), where 

G/.(e) = (G/,(e)"i-f E/.(e))-i. (9) 
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Figure 1: (a) Comparison of the f-electron density 
of states calculated by the exact diagonalization 
of the finite system (full lines) and by IPX (broken 
lines) for J/ = 0, 1, 2, 3 (from bottom to top), 
(b) The f-electron density of states calculated by 
SCSOPT for the same parameters as in (a). Note 
that the origins of the vertical axis are shifted. 

3 Results 

It is rather easy to solve the above equations self- 
consistently when h — 0. We choose W = 1 
as a unit of energy and V — 0.5. The fol- 
lowing calculations are limited to the absolute 
zero temperature. In Fig. 1(a), we compare the 
f-electron density of states p(e) by IPT for [/ = 2 
with that obtained in I by EXD using the finite 
system, There are rich structures in p{e) ob- 
tained by EXD. In addition to the hybridization- 
split Kondo peaks at small energies, there are 
peaks at e ~ ±0.8 and ±1.7 for U = 2, for in- 
stance. The spectra by SCSOPT are displayed 



Figure 2: The real and the imaginary parts of 
the f-electron self-energy calculated by IPT (full 
and broken lines for real and imaginary parts) 
and SCSOPT (dash-dotted and dotted lines) for 
U=2. 



in Fig. 1(b). The whole structures are better re- 
produced by IPT than SCSOPT. This is in ac- 
cord with the results obtained by Mutou and 
Hirashima|22 on the comparison of the quantum 
Monte Carlo calculation with IPT and SCSOPT. 
The peaks at e ~ ±1.7 in EXD are the lower and 
the upper Hubbard bands, which are pushed out 
a little beyond ±U/2. In SCSOPT, there are only 
small shoulders at these positions, but it is inter- 
esting that there are larger peaks at e '--^ ±0.8 
which correspond to the peaks in EXD at the 
same positions. 

These differences in the spectra in SCSOPT 
and IPT come from the differences in the self- 
energies. We have plotted the f-electron self- 
energies for [/ = 2 in Fig. 2. It is seen that the ab- 
solute values of the real and the imaginary parts 
of I]/(e) are larger in ITP than in SCSOPT. This 
makes the Hubbard bands in IPT be positioned 
at the right energies and enhanced sufficiently to 
reproduce the EXD results. The fine structures 
of the spectra in SCSOPT seen above seem to 
be originating from the structures in the imagi- 
nary part of the self-energy at the corresponding 
energy regions. 

Despite the better description of the f-electron 
density of states in IPT than in SCSOPT, we 
have found that magnetic properties are better 
described by SCSOPT (see below). This is be- 
cause IPT is not a conserving approximation. |^| 




Figure 3: (a) m vs. mo curves (see the text) 
obtained by SCSOPT for [/ = 2 and h =0, 0.1, 
0.2, 0.3, 0.4, 0.5 from bottom to top. (The IPT 
result for ft, = is displayed by the dotted curve 
for comparison.) (b) The total energies for the 
same set of parameters are plotted from top to 
bottom. 



Therefore, we will use SCSOPT in the following 
calculations. 

Under finite magnetic field, one also has to de- 
termine n/o- self-consistently. We found, however, 
that the convergence of the iterative calculation 
becomes rather poor in the vicinity of the criti- 
cal field. To overcome this difficulty, we first as- 
sume appropriate input values n'j^ for Ufa- and 
solve the above equations for fixed Ef^. Then 
we obtain Gf„{e) and nf„ from them. We re- 
peat this calculation for various rflr^s. and plot 
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in Fig. 3(a). 



tion 771 = mo. It is clear from the S-like shape 
of m vs. mo curves that the first order transition 
takes place at a certain critical field he, which was 
found to be he — 0.101 &iU — 2. We also checked 
that the total energy, p3[ 



Etc 
where 



j defie) ( - ) Im eGe.(e) 
1 



e--^f,ie)]Gf4e) , (10) 



ek-Vyie-Ef,-^f4e)) 



F{e-V^g,ie)), 



(11) 



The points of intersection with the line of the 
slope of 45 degrees give the self-consistent solu- 



is lower for larger mo (Fig. 3(b)). 

The transition becomes continuous for U < 1.7 
since the S-shapc of the m vs. mo curves becomes 
weak. On the other hand, the ferromagnetic state 
becomes stable even at ft = for U > 2.6 because 
of a stronger S-shape. 

The data in Fig.|^ are spline-interpolated on 
mo and ft, from which we determine the self- 
consistent solution m = mo with a good accu- 
racy. We plot the magnetization curve in Fig. 4. 
It shows a jump at he for 1.7 < U < 2.6, although 
the size of the jump is larger than that in EXD. 
The double steps found by EXD ioi U — 1 could 
be a reminiscence of a change of the character of 
the transition from first to second order. 

IPT also yields S-shape curves in the m vs. 
mo relations, but the S-shape is stronger than in 
SCSOPT (see Fig.3(a), the dotted curve), giving 
smaller he {he — 0.071 for U = 2). A fatal draw- 
back in IPT is that it gives wrong value for the 
initial slope of the magnetization curve (i.e. mag- 
netic susceptibility Xf)- We found Xf =1.22 and 
1.20 in SCSOPT and EXD, respectively, whereas 
Xf = 0.67 in IPT for U = 2. This discrepancy 
is due to the fact that IPT is not a conserving 
approximation ||l^ as mentioned above. 

The f-electron densities of states under finite 
magnetic field are plotted for [/ = 2 in Fig. 5, 
from which the energy gaps are obtained and are 
plotted in Fig. 6 as a function of the field. The 
gap closes abruptly between ft =0.101 and 0.1015. 
From these results, we consider that the first or- 
der transition obtained in I may not be due to the 
finite-size effect. 

Finally, we would like to comment on the be- 
havior of the self-energy in the magnetic field. 
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Figure 4: Magnetization curves for U =0 (dot- 
ted line), 1 (triangles), 2 (squares), 3 (diamonds) 
calculated by SCSOPT (filled symbols) and EXD 
(open symbols). 





Figure 6: The Energy gap obtained by SCSOPT 
(circles) as a function of ft- for f7 = 2. The dotted 
line indicates the gap for U = scaled a.t h — 0. 



Fig. 7 shows the field dependence of the real part 
of the self-energy. As is seen clearly, KeT,fa{e) 
does not depend on h when h < he- This is be- 
cause the shift of Paif-Ys in eq.(||) by the magnetic 
field does not yield a noticeable change of ^/^(e) 
until the gap is closed by the sufficient field. Due 
to this behavior in the self-energy, the renormal- 
ization factor z changes little for h < he- Thus the 
gap closing is mainly caused by the relative Zee- 
man shift of up- and down-spin bands. Note that 
the gap does not close like Eg{h) = Eg{Q) — 2h 
even for [/ = since we apply magnetic field only 
to f-electrons, so that the shape of Pcr(e) canges 
with the field. The energy gap for ft, = is 
given by Eg{0) = -W + VW^TW^ = 0.414 
for W — 1 and V — 0.5, whereas ftc is given by 
he = V^/W = 0.25, which is larger than Eg{0)/2. 
We have plotted Eg{h) for /7 = in Fig.6 by the 
dotted line, which is scaled at ft = 0. The gap 
by SCSOPT decreases a little faster than that for 
U ~ 0, which means that the shape of Pai^) at 
U = 2 changes more than that of [/ = by the 
field despite the small change of z. 



Figure 5: The f-electron density of states for up- 
spin electrons {U = 2) under the magnetic field 
ft =0, 0.05, 0.101 and 0.1015 from bottom to top. 



4 Discussions 

We have found that the energy gap in Kondo insu- 
lators is closed by the application of a sufficiently 
strong magnetic field of the order of the gap, and 
a transition from insulator to metal does take 
place in the limit of infinite spatial dimensions. 
This is in strong contrast to the suggestions for 
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Figure 7: The real parts of the self-energy are 
plotted for h =0 (full line), 0.05 (dashed line), 
0.101 (broken line) and 0.1015 (dotted line). 



the one dimensional Kondo insulators. |jT^ |l^ We 
found that the transition is of first order, and the 
magnetization shows a jump at the critical field 
/ic, irrespective of the methods of calculations, al- 
though the first order transition is limited to the 
range 1.7 <U < 2.6 in SCSOPT. The calculation 
was only for the electron-hole symmetric case. It 
is of much interest to extend the present calcula- 
tion to more general situations and to more real- 
istic models with orbital degeneracy. We started 
from a paramagnetic ground state, assuming that 
a frustration should have destroyed a long range 
order. 1 21 It is necessary, however, to extend the 
calculation to such an explicit model which has a 
paramagnetic ground state. 

Despite these issues to be clarified in future 
studies, we stress that the magnetic-field-driven 
insulator-to-metal transition is of great interest 
and importance for the study of the characters 
of the Kondo insulators since a sufficiently strong 
magnetic field always tends to close the gap. On 
the other hand, the pressure, for instance, some- 
times increases the mixing and the gap, but in 
some other cases increases the band width, caus- 
ing the overlapping of the band and rendering 
the system semimetallic. In this sense, the mag- 
netic field is much simpler and a good tool for 
the study of the insulator-to-metal transition in 
strongly correlated electron systems. 

First order transition has not yet been observed 
in YbBi2. This is partly because the energy 
gap was determined from the temperature depen- 



dence of the resistivity, so that the finite tempera- 
ture might have smeared out the transition, or the 
system is out of range of the first order transition. 
More direct experiments at low temperatures, e.g. 
measurement of the dynamical conductivity, etc. 
are desirable. Lacerda applied the magnetic field 
of up to 50T to SmBg and did not find any indi- 
cation of a transition to metal. |Q This may be 
because of the smaller magnetic moment for Sm 
{gj — 2.14 for Ce whereas gJ = 0.71 for Sm). 
Magnetoresistance of CeNiSn is measured to be 
positive]^ in contrast to the present calculation. 
The structure of the gap and the mechanism of 
the low temperature transport in this compound 
is still not clear. However, the effect of the mag- 
netic field on the strongly correlated state should 
be different if the gap is of V shape. Study of such 
situation will be interesting also from the present 
theoretical point of view. 

Recently, a new IPT scheme was proposed [p5|, 
^ , which can be applied even to the case without 
electron- hole symmetry. It may be applicable to a 
calculation of the properties of the strongly corre- 
lated systems under magnetic field, and may im- 
prove the present calculation, including the cases 
with larger U and the study of metamagnetism 
in the metallic heavy fermion systems. 
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